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Abstract 

This paper proposes a novel joint non-binary network-channel code for the Time-Division 
Decode-and-Forward Multiple Access Relay Channel (TD-DF-MARC), where the relay linearly 



^ combines - over a non-binary finite field - the coded sequences from the source nodes. A method 

1^ based on an EXIT chart analysis is derived for selecting the best coefficients of the linear combi- 

nation. Moreover, it is shown that for different setups of the system, different coefficients should 
pg be chosen in order to improve the performance. This conclusion contrasts with previous works 

L| where a random selection was considered. Monte Carlo simulations show that the proposed scheme 

/\ outperforms, in terms of its gap to the outage probabilities, the previously published joint network- 



channel coding approaches. Besides, this gain is achieved by using very short-length codewords, 
which makes the scheme particularly attractive for low-latency applications. 
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I. Introduction 

Node cooperation has been widely shown to improve the performance of wireless networks 
with several terminals by increasing the robustness of the system to channel variations such as 
deep fades, as well as by enabling significant energy savings. The essence of node cooperation 
lies in jointly processing the in-network information by all the constituent nodes of the 
network, which allows improving the spectral and power efficiency of wireless networks and 
ultimately, attaining the desired diversity-multiplexing tradeoff without requiring additional 
complexity (e.g. co-located multiterminal MIMO schemes). 

A particular example of relay cooperation in multi-terminal networks is the so-called 
Multiple- Access Relay Channel (MARC). The MARC is a communication scenario where 
two or more information sources forward data to a single common destination with the help 
of an intermediate relay [|T||. In this scenario, the relay can work in full [[T|, [|2| or half duplex 
mode. For the half duplex mode, the following transmission strategies have been proposed: i) 
the Constrained MARC (C-MARC [[3|), where the sources transmit during the first time slot 
and coordinate with the relay during the second time slot by transmitting information; ii) the 
Orthogonal MARC (0-MARC, see Q, [|5|), where the sources and the relay transmit over 
two orthogonal channels; and iii) the Time-Division MARC (TD-MARC, see [6|), where 
both sources and the relay convey their data by using three orthogonal channels, i.e. for two 
sources the total transmission time is divided into three time slots, one for each transmitting 
node. On the other hand, during the last decade the well-known relaying strategies Decode- 
and-Forward (DF), Compress-and Forward (CF) and Amplify-and-Forward (AF) originally 
developed for the conventional relay channel [|7| have been applied to the aforementioned 
MARC models in a number of contributions (see e.g. |[l|~l|l| ^nd references therein). 
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Fig. 1. Block diagram of the considered 2-user TD-DF-MARC scenario. 



When dealing with practical coding schemes for the MARC scenario, the TD-MARC, 
along with a DF strategy, has been widely studied in the related literature ||6|, [[8|-p7|. 
Although the time-division scheme involves a suboptimal use of the available bandwidth, it 
allows for an easier implementation in practical systems thanks to the use of half-duplex 
relays and the lack of stringent synchronization constraints. Besides, the DF strategy offers 
a higher code design flexibility. In this work we specifically focus on the 2-user TD-MARC 
with a DF relaying strategy (hereafter coined as TD-DF-MARC), which is schematically 
depicted in Figure [T| It is important to note that the capacity bounds for this model can be 
derived from the capacity bounds of the 0-MARC p|, p^|. 

In the MARC we are interested in maximizing the information conveyed by the relay node. 



To this end, network coding [19|, [20 1 has become a widely used technique to complement 
channel coding schemes used for combating channel-induced errors. Combining the data from 
both sources at the intermediate relay node embodies a practical tool for approaching the 
capacity bounds of the TD-DF-MARC scenario. However, by treating network and channel 
coding separately some performance loss is expected, since the network decoder cannot use 
the output soft information computed by the channel decodes. Likewise, the channel code 
cannot exploit the redundant information provided by the network code. This observation is 



further supported by the findings in [21 1, [22[, where it was shown that in general, capacity 



can only be achieved by conceiving channel and network coding as a single non- separated 
data processing stage. 

In this context, several practical joint network-channel coding schemes have been recently 
proposed [[6|, [[8|-[[17|, ||23|. In [[6| joint network-channel coding for TD-DF-MARC model 



was first considered. The authors proposed distributed regular LDPC codes as the joint 
network-channel code at the relay node, where the destination jointly decodes the messages 



from the sources with the aid of the information sent from the relay, as opposed to p4 | 
where the two messages transmitted from the sources were separately decoded. In [[8| the 
authors follow the same idea by proposing a turbo-code-based joint network-channel coding 
scheme. Parallel to these proposals, the authors in [9] proposed a similar scheme for high- 
order modulations. More recently a joint coding scheme based on WiMax LDPC codes was 



presented in [10], whereas in [11|, [12| two schemes based on turbo codes were investigated. 



Furthermore, some joint non-binary coding schemes have been recently reported in p3| (non 



binary network coding) and [16|, [17| (non-binary network and channel coding). Finally, in 



1 13 1, [14 1 the authors proposed a joint network-channel coding scheme where the relay 
transmits the soft values resulting from its decoding procedure over AWGN channels. 

The present work joins the upsurge of research on the TD-DF-MARC scenario by proposing 
a novel Joint Network-Channel Code (JNCC) where the relay linearly combines - over a 
non-binary finite field - the coded sequences from the source nodes. The iterative decoding 
procedure at the common destination is performed by running the Sum-Product Algorithm 



(SPA [25 1) on the factor graph describing the proposed JNCC, which is compounded by three 
sub-factor graphs: two describing the channel codes of each source, and a third describing 
the network coding operation performed at the relay node. Specifically, the key contributions 
of this manuscript over the state of the art on this topic are as follows: 

• The proposed scheme does not perform channel coding on the already network-coded 



bits, reducing the complexity at the relay node without compromising performance. To 
the knowledge of the authors, all practical schemes for fading channels found in the 
related literature^ perform channel coding on the already network-coded bits. 
It is shown that a tailored selection of the set of coefficients used in the network coding 
operation, namely, the Network Coding (NC) coefficients, outperforms a random choice 
as done in (T6\. This selection is performed by matching the EXtrinsic Information 



Transfer (EXIT, [26|, [27|) functions of the compounding codes on the EXIT chart. 

• Contrary to previous hterature, the JNCC allows the sources to use completely different 
channel codes, at the sole expense of an increased complexity when choosing the NC- 
coefficients through EXIT-curve matching. 

• The work presented here considers convolutional codes at both source nodes. As a conse- 
quence of the previous point, both convolutional codes can be independently terminated. 
This fact allows us to use very short-length codewords, making the scheme particularly 
attractive for low-latency applications. 

The remainder of the manuscript is organized as follows: Section [ll] introduces the system 



model, whereas the decoding algorithm of the proposed JNCC is detailed in Section III In 



Section IV an analysis on the influence of the NC-coefficients is performed through EXIT 
charts. Section |V] discusses the obtained Monte Carlo simulation results, and finally Section 
l^ends the paper by drawing some concluding remarks. 

II. System Model 

Let (n, (3, V) be the underlying probability space where all the random variables (r.v.) are 
defined. We use uppercase when referring to r.v. and lowercase when referring to realizations 

'in [131 a similar coding procedure at the relay is proposed; however, their study is restricted to AWGN channels and 
binary network coding. 



of r.v. In addition, we use boldface when referring to vectors; thus, uppercase and boldface 
refer to random vectors. For discrete r.v., we denote the probability mass function (p.m.f.) 
of the discrete r.v. X as Px{x) = V{X = x}. For continuous r.v., we denote the probability 
density function (p.d.f.) of the continuous r.v. X as pxix). However, when the context is 
clear, we use P{x) and p{x) for p.m.f. and p.d.f., respectively. 

Referring to Figure [TJ for simplicity we have considered a symmetric scenario consisting 
of 2 unit-entropy binary information sources Si and ^2, which generate blocks U^ G {0, 1}^ 
and U^ G {0, 1}^ of length K. As depicted in this figure, at each transmitter the sequence is 
channel-coded by a convolutional code, producing the codeword C" = {C™}^]^ G {0, 1}^, 
with m E {1,2} denoting the source index. The code rate is therefore given hy R = K/N . 
Each codeword is then interleaved yielding the interleaved codeword X™ = Ilm (C™), where 
Hi and 112 are two different random spread interleavers with a spread factor equal to g G N. 
Finally, the codeword is modulated, resulting in the transmitted sequence S™ = {Sl"'}f^i, 
which is transmitted over M = N/2 complex dimensions (i.e. A^ real dimensions). During 
the first and second time slots source Si and 52 transmit to both the relay and destination 
nodes the sequence S^ and S^, respectively. The third time slot is used by the relay to process 
the data from the sources and transmit the resulting coded data to destination. 

Regarding the links between nodes, we denote as dkj the distance from transmitter k G 
{1,2, R} (R: Relay) to receiver j G {R, D} (D: Destination). Moreover, considering the 
power at the end of the source-destination link Pq as the reference, the received power at the 
end of each link will be given by Pq ■ {ds,D/dk,jY , where 5 denotes an attenuation exponent. 
In what follows, and without loss of generality, the distances are normalized with respect to 
ds^D = 1 and we consider Pq = I. Thus, the attenuation undergone by the signals due to the 
distance-dependant propagation losses of a given link can be expressed as (i^l Therefore the 



received symbol per real dimension at each receiver is given by 

Yr = a''^ ■ ^^^ ■ S^ + N':'\ (1) 

where a^j is Rayleigh distributed with -E'[afcj] = 1 V/c, j, and {N^ '■^}^i are modelled as real 
Gaussian i.i.d. random variables with zero mean and variance N^'-' /2. The values of {0;^.,} 
are assumed to remain constant within the duration of a transmitted block (i.e. quasi-static 
fading). Moreover, full channel state information (CSI) is assumed at the receivers. 

A. Relay Node 

Consider the set of all 2'' polynomials p{z) of degree g — 1 with coefficients lying in GF{2) 
(the binary Galois field). Let g(z) be a prime polynomial (i.e., monic and irreducible polyno- 
mial) of order q. Then, this set becomes a finite field, GF(2'^), by defining the addition © and 
multiplication (g) rules as the mod g{z) remainder of the sum and product of two polynomials, 
respectively. Notice that since the mod g{z) addition rule is just componentwise addition of 
coefficients in GF(2), GF(2fl) under addition is isomorphic to the vector space {GF(2)y of 
binary g-tuples with mod 2 elementwise addition, denoted hereafter as A. Therefore, there is 
a one-to-one mapping il)q : {GF{2)y — > GF{2'^) defined as ^jJq{ao, . . . ,ag_i) = J2k=o'^kZ^ 
such that i{jq{ai) © V^g(b) = ipq{sL A b), where a, b G {GF{2)y. In addition, we index the 
elements pi E GF{2'^), i G {0, . . . , 2"^ — 1} by the base-10 notation of the corresponding binary 
tuple (oo, . . . , ttq-i). In the following we refer to the elements of the finite field GF{2'') as 
non-binary symbols. 

In the first and second time slots the relay receives the channel sequences Y^'^ = {Y^ ' }^^ 
and Y^'^ from sources S^ and S"^, respectively and it deinterleaves them. Then, it executes the 
BCJR algorithm [28 1 twice in order to obtain the estimations C^ and C^ of the source channel- 
coded sequences, which are then interleaved in order to obtain the estimated interleaved coded 



bits Xi and X^. 

Each of the interleaved coded sequences X^ and X^ is partitioned into N/q sub-sequences 
of length q. We denote as V^ = V'q({-^l"}i=i) ^ GF{2'^) to the non-binary symbol associated 
to the corresponding sub-sequence. The non-binary symbol of the relay V^ is now computed 
as the linear combination of the non-binary symbols corresponding to each source, i.e. 

where h = {h^,h'^) and /i™ G {pi}"^!^^ represents the NC-coefficients used in the linear 
combination. Finally, the modulated symbols associated to each sub-sequence are computed 
as {Sf^}f^^ = 2 ■ i>~^{V^) — 1 and the transmitted signal S^ is obtained by concatenating 
the N/q resulting modulated sub-sequences. 

III. Iterative Joint Network-Channel Decoder 

The destination receives the channel outputs Y = (Y^'^, Y^'^, Y^'^). The aim of the 
JNCC decoder is to find the source binary symbols {f/™}^!^ that maximize the conditional 
probability P(M™|y), which is obtained by marginalizing the joint conditional probability 
P(u™|y). This marginalization is efficiently computed by applying the SPA over the factor 
graph describing P(u'^|y). Figure [2] shows the three compounding sub-factor graphs of the 
proposed JNCC: two describing the source convolutional codes, and the third one describing 



the network code used at the relay. As explained in Section II-A the factor graph of the 
relay network code is in turn composed of N/q parallel and identical factor nodes, depicted 
in Figure |2] as the oversized factor nodes which we hereafter refer to as network check nodes 
and labelled with N'Ci, / G {1, . . . , N/q}. Furthermore, we define y/ — (y^ ' , y^ ' , y^ ' ) as 
those components of Y associated to the network check node AfCi. In the next subsection 
the derivation of the factor graph corresponding to one of these network check nodes is 



explained, which is then incorporated into the overall factor graph plotted in Figure [2j 



Convolutional Code [S^] 



Fig. 2. Factor Graph of the Proposed JNCC. 
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Since the overall factor graph of the JNCC has loops, the SPA is iteratively run between 
the sub-factor graphs corresponding to the relay network code and the convolutional codes. 
After a fixed number of iterations X, the f/™ is computed as 



P«|y) (X Y, Tkisk, u^, C-, s,+,)a{s,)P{sk+i) J] ^(^D, 



(3) 



where a and /3 are the forward and backward messages passed from the adjacent state nodes 
to the factor node Tk given by the Trellis of the convolutional code; c™ are the coded bits 
c™ associated to n™; and 7 are the messages passed from the variable nodes c™ to T^ (i.e. 
the likelihoods). Note that in the case of not having a relay, the likelihoods 7 are given by 
7(c™) = p{y^' Ic^)- However, when the relay is present, these likelihoods now depend on 
the messages passed by the network check nodes associated to the interleaved binary symbol 
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x"?.!, ,, i.e. 7(c™) = 7(x"?_i, A where 7(x™) oc p{yi\x'^). In the next subsection a factorized 
form of p(y;|x™) is derived. 

A. Factorization of p{yi\x'^) 

In this section the factorization of the conditional probability piyi\x]^) is derived. As 
before, the upper index j refers to the random variables associated with both sources and the 
relay, whereas the index m only refers to variables associated to the sources, i.e. j E {I, 2, R} 
and m E {1,2}. Besides we use index m to refer to the complement of m, i.e. m = 3 — m. As 
p{yi\x]^) only depends on the bits belonging to the corresponding pair of sub-codewords of 
length q, we focus on any given pair {X"^}m=i 2, where X™ = (X[", . . . , Xjp). Thus, we use 
the subindex i to denote the position of a given bit inside its corresponding sub-codeword. 
For the sake of simplicity we drop the super index D from the set of received signals 
Y = {Y^'^}j=i 2,ij and rewrite the channel random variables for a given sub-codeword as 

In order to compute piyi\x'1^), let us first focus on the a priori information provided by 
the channel decoder. The sequence of random variables {Xp}j^^ is assumed to be i.i.d. 
based on the fact that a spread interleaver is used to suppress strong dependencies (4-length 
cycles in the underlying factor graph) between the bits belonging to the same sub-codeword. 
Moreover, the p.m.f of its associated non-binary symbol V"^ = tpg{X^, . . . , X™) is given by 

PvAv)= E ^[v = U<,---,^7)]IlPxri^T), (4) 

a;5",...,a;™ i=l 

where the last factor represents the a priori probabilities of {X^}f^^. 

As shown in expression (|2]), the non-binary symbols V^ and V'^ are linearly combined 
over the finite field GF{2''), producing the non-binary symbol V^. Therefore, the non-binary 
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symbols joint p.m.f. can be factorized as 

P{v\v^,v'^) = P{v'^\v\v^)P{v^\v^)P{v') = P{v')P{v^)l [v^^ = f{v\v^)] , (5) 

where V^ = f{V^,V^) is the linear combination defined in (|2]). Now let us focus on the 
information coming from the channels. Due to the TDMA scheme, the following Markov 
chains are verified, Y^ ^ V^ ^ (F^ 1/^, Yf, Yf ); Yf ^ V^ ^ (F\ V^, Y/, Yf ); 
Yf ^V^^ {V\ V\ Y}, Yf). Hence, we have 

p{yi\v\v^v'') = p{yl\v')p{y^\v^)p{yf\v''). (6) 

Furthermore, the non-binary symbols can also be expressed by the corresponding modu- 
lated symbols {SljU e {il}"^ as V^ = ^, ({(1 + Si)/2}Q ^ 9, {{S^U), yielding 

piylW) oc Y. 1 [^' = ^<?(^i' • • • ' 4)] rip(^{>'), (V) 



•^1 V -i^q 



where the last product is due to the memoryless channel assumption made in this work, with 



Piyii\4) ocexp 






(8) 



Now, from expressions (|5]) and ([6]) and by applying the Bayes theorem, we obtain the joint 
a posteriori p.d.f. of the non-binary symbols as 

p{v\v\v^\y,)o,p{y]\v')p{y'Av')p{yf\v^) ■![.« = f{v\v')] ^"^^J^^^^j'^^'^ (9) 

Likewise, one can compute the a posteriori p.m.f. P{v'^\yi) of the non-binary symbol 
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associated to a given source by marginalizing the previous equation. Thus, 

On the other hand, following a similar reasoning as in ([4]), we obtain 

from which we can compute the bitwise channel conditional p.d.f. by marginalizing and 
applying again the Bayes Theorem, i.e. 

p{yi\^T)= E ^(^™ly0iK = V^.«'---'O]|7f^' (12) 



where ~ x™ = {xf}\jj^i. Finally, combining ([10]) and ([12]), we get 



p{yi\xT)= E i[^' = ^,«,---,0]-p(yrK)n^"(^")^""(^")' (^3) 



with 



pu..c^^m^ A y- 1 [^/? ^ /(t;\t;^)] ■p(yr>"^)P"(^^"^).p(yf b""). (14) 



The factorized form of p{yi\x'^) given in expression ( [13] ) is graphically represented by the 
factor graph depicted in Figure [3} where for the sake of clarity the sub-index / is dropped. The 
application of the SPA over the factor graph of Figure [3] allows for a efficient computation 
of the likelihoods 7(x™). To be concise, the SPA iterates between the sub-factor graphs 
corresponding to the relay network code and the convolutional codes. It should be remarked 
that if the probability p(yf |f ^) (message) is not dependent on the data from source m (e.g., 
when the NC coefficients are set to zero or the relay-destination channel is in deep fade), 
then pMARc^^m-j ^in i^g uniformly distributed and consequently the exchange of messages 
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between both convolutional decoders (p^ARcj^^m-j^ ^^^ ^^^ improve the performance of the 
decoder (see Fig. [3]). 

Although the factor graph shown in Fig. [2] has been constructed for convolutional codes, 
it could be easily modified if iteratively decodable codes (e.g. LDPC, Turbo) are used, by 
facing their outer codes with the Relay Network code subgraph. However, care should now 
be taken when programming the decoder activation scheduling. 



In [29 1 scheduling algorithms for both parallel and serially concatenated codes with several 



compounding graphs are proposed. We further refer to [29| for algorithms that find the fastest 
convergent code activation schedule. Finally, a particular case where non-binary LDPC codes 



were used as channel codes was published by current authors in [17|. 
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Fig. 3. Factor Graph representing the factorization of p{yi\xl"), i.e. the network check node AfCi. 
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In the next section, the EXIT charts of the proposed JNCC are used to obtain a pair of 
NC-coefficients that optimizes the performance of the system. To simplify the exposition, the 
optimization procedure next presented assumes that both sources are detected with no errors 
at the relay, in line with the assumptions made in the related literature [|6|, [|9|, [10|, [17|. 



IV. Analysis of the Relay Network Code based on EXIT charts 

An additional insight can be gained through an analysis of the NC-coefficient based on the 
interchange of mutual informations between the output of the channel decoders and the output 
of the relay network decoder. Due to the iterative nature of the decoding algorithm, EXtrinsic 



Information Transfer (EXIT, see [27|) charts are a good method for visually exploring this 
iterative exchange of information. Given a code, the EXIT function associated is defined by 
the relation between the a priori mutual information at the input of the decoder (commonly 
denoted as la) and the corresponding extrinsic mutual information I^, at its output, i.e. /g = 



T{Ia). For further information on EXIT charts we refer [27 1 and references therein to the 
reader. 

We will denote the transfer function of the network code for a given q and h as Jg"^ = 
^h(-^D- Notice that, T^(0) and T^(l) represent the extrinsic information at the output of the 
network decoder with no a priori and full a priori information about the information bits, 
respectively. Moreover, since the mutual information at the input of the channel decoders 
is equal to the mutual information at the output of the network decoder, i.e. I'^ = I^'^, 
the extrinsic mutual information at the output of the convolutional decoders is given by 
jcc _ j^ccj^jNc-j xhus, for a successful decoding procedure, there must be an open gap 
between both EXIT curves so that the iterative decoding can proceed from I^'^ = to 
Jg*^ = 1q When both transfer functions intersect, the iterative process will stop at a given 

^Although for a perfectly successful decoding the final J^'' should be equal to one, we also consider the values of I^'^ ~ 1 
that yield to negligible error floors. 
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extrinsic mutual information of the source bits Jg*^ < 1. A crossing yielding Jg*^ < 0.5 will 
be referred to as early-crossing and as late-crossing, otherwise (i.e., 0.5 < Jg'^ < 1). Since 
the transfer functions are monotonically increasing, the higher the value T^(0) is, the later 
the early-crossing will occur. On the other hand, and if no early-crossing occurs, the higher 
the value of T^(l) is, the closer I^"^ will be to one. 

Next, we justify the reason why the NC coefficients have a greater impact on the value of 
T^(l) than on T^(0). When no a priori information is available at the network decoder, all 
its information comes from the sources- and relay-destination links. On the other hand, when 
full a priori information is available, only the information provided by the relay-destination 
link is relevant, since the information regarding the coded messages from the sources is 
fully supplied by the channel decoders. In addition, from the Area Theorem of EXIT charts 



[27 1, which states that the area below the transfer function depends only on the rate of 
the encoder, the area below the transfer functions will be constant Vg, h. Consequently, the 
transfer functions of those NC-coefficients with a large value of T^(l) are expected to be 
flat shaped at low values of 1^"^, and steep shaped at values near one. 

To corroborate the above, let us first assume AWGN channels (i.e. ol = (1, 1, 1)). Figure |4] 
plots the EXIT chart of a network check node for different NC-coefficients and values of q, 
along with the transfer function of a 2^-state (CC2) and a 2^-state (CC6) convolutional code 
with same transmission rate 1/3. Note that by the symmetry of the network, it is sufficient to 
consider only those NC-coefficients h = (pi, Pj) with indexes verifying l<i<j<2'^ — I. 
For the sake of clarity, we will denote h = {pi, Pj) by {i,j). From this figure, it can be 
observed that the transfer function corresponding to h = (1, 1) (i.e. raw XOR-network coding) 
maximizes T^(0) and minimizes T^(l) for all q. Also observe that for high values of q and 
for some NC coefficients, the initial flat shape of the plotted transfer functions curves induces 
an early-crossing with the curve of CC2, and consequently the system requires more energy 
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Fig. 4. EXIT charts for different network decoders and convolutional codes. 



per symbol to open the gap between these curves. In this case, one may opt to select NC- 
coefficients that increase the value of T^(0), so the early-crossing could be avoided. As a 
drawback, a higher error floor due to their associated lower T^(l) is expected. As highlighted 
in the figures, such alternate coefficients are given by h = (2, 2) {q = 2), h = (6, 6) {q = 3), 
h = (7, 7) (g = 4) and h = (14, 15) (g = 5). 

Therefore, for AWGN channels we can conclude that by a proper choice of the coefficients 
an iterative gain will be achieved when using a linear combination of the estimated symbols. 



as it was also stated in [ 13 1. This contrasts with the more commonly used scheme based on a 
channel encoder at the relay. Moreover, from the variety of transfer functions generated by the 
family of NC coefficients, it is also concluded that a tailored non-binary linear combination 
at the relay can outperform 1) the XOR coding method (i.e. h = (1, 1)) first proposed in 



1 11 1, [13 1, [24 1; and 2) a random choice of such parameters as was proposed in [15|, [16|. 
Before considering the more general case of fading channels, observe that the proposed 
soft-output network decoder can be seen as N /q parallel systematic codes W(h); each of them 



17 



having the 3g-length codewords W = {{XDI^^, {Xf}1^^, {-^j^lLi) where its systematic 
bits, ({X/}^^^, {Xf}1^^), correspond to the interleaved bits of the sub-sequences from sources 
S^ and S"^, respectively, and its parity bits, {X^}'^^^, correspond to the bits generated by the 
relay. Hence, each pair of NC-coefficients generates a particular code W(h) with a different 
distance spectrum and therefore, different transfer functions. Moreover, the relation among 
these functions depends solely on their distance spectrum and it is independent of the quality 
of the channels. 

We now look at Rayleigh fading links. In this case, different realization of the coefficients 
of the Rayleigh fading a will produce differently shaped transfer functions, making T^ a 
random mapping from M — ^ M. If the channel between the relay and the destination is in a 
deep fade, no gain is obtained by iterating (since pi^ARc^.^ j^ uniform), and a late-crossing 
might occur (given that T^(0) is high enough). On the other hand, if a deep fade occurs in 
both source-destination channels, an early-crossing could be produced due to the low value 
of T^{0). Both situations are shown in Figure Bl which plots the EXIT curves of 100 channel 
realizations for g = 2 and h = (2, 3). In this Figure, the curves corresponding to two fading 
realizations producing an early- and late-crossings, have been highlighted for clarity. Since a 
decoding error is produced when a crossing occurs, the probability of a failed decoding event 
could be approximated by Pjfailed decoding} f^ "Pjearly-crossing} + P{late-crossing}. 

It should be mentioned that the Pjearly-crossing} strongly depends on the quality of the 
channels whereas P{late-crossing} is not so dependant. The reason being that increasing 
the received signal-to-noise ratio will reduce the probability of an early-crossing, since 
the decoder would be fed during the first iteration with more reliable channel information. 
However, the late-crossing probability will not be significantly reduced since the influence 
of the information provided by both source-destination channels decreases as the channel 
decoders begin to provide a priori information regarding the source encoded bits. Therefore, 
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as the signal-to-noise ratio increases, the late-crossing probability starts to dominate the failure 
probability, and hence, beyond a certain value of the signal-to-noise ratio, the early-crossing 
probability becomes negligible, regardless of the channel code used. 




jNC jCC 



Fig. 5. EXIT curves of 100 system usages for q = 2 and h = (2, 3). 



In addition, the failure probability will also depend on the type of convolutional code 
used and as shown in Fig. Isl the more complex the code is, the larger the value of T^(0) 
should be in order to avoid early-crossings. On the contrary, the less complex the code is, 
the higher the value of T^(l) should be in order to avoid late-crossings. Therefore, the value 
of T^(0) and T^(l) will strongly determine the failure probability of the proposed decoder. 
To analyze the crossing probabilities one has to statistically characterize the behavior of the 
random variables T^(0) and T^(l) for different values of q and NC-coefficients. To that end, 
we next show that the distribution of these random variables depend on the path-loss gain 
suffered by the signal coming from the relay (i.e. the position of the relay). 

For T^(l) (i.e., the full a priori case) the extrinsic information generated by the network 
check nodes (i.e., P^'^'*'^(-)) depends on the quality of the relay-destination channel through 
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the check node associated to the linear combination. Therefore, in the limit when the SNR of 
the relay-destination channel tends to infinity, the probability of V(T^{1) = 1) tends to one. 
As a consequence, the slope of transfer functions will increase since T^(0) remains fairly 
constant; and consequently, some of the early-crossings will be avoided. 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 6. CDF of T^(0) and r«(l) for g = 3 and all possible values of h. 



To make this analysis clearer. Figure [6] shows the Cumulative Distribution Function (CDF) 
of T^(0) and T^(l) for g = 3 and all combinations of h for the cases where: the relay and 
the sources are at the same distance to the destination (0 dB gain), and the relay is placed at 
half distance between the sources and destination (10 dB gain). As mentioned for the AWGN 
case, it can be observed that the choice of a different value of h has a stronger impact on 
T^(l) than on (T^(0). Also, increasing the quality of the relay-destination channel leads to an 
increase in the values of the realizations of T^(l), which in turn diminishes V{T^{1) < x} 
for any x E [0, 1]. Moreover, V{T^{0) < x} does not significantly change as the gain in the 
relay-destination link is varied. 

Observe from Figure [6] that given a q, there is a particular value h^ (shown as solid lines) 
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that minimizesj'P{T^(l) < x}, Vx G [0, 1] regardless the quality of the channels. Similarly, 
there is a value h^ (shown as dashed lines) that minimizes V{T^{0) < x}, Vx G [0, 1]. 

At this point, by denoting the smaller SNR value that makes the term P{late-crossing} to 
dominate the system performance as SNRth, the best choice of NC-coefficients is h^ when 
SNR > SNRth, and hg, otherwise. Based on this, the steps for selecting the NC-coefficients 
can be summarized as follows: 1) find the values of hg and h*; 2) estimate by Monte Carlo 
simulations the value of SNRth', and 3) choose Iiq or h^ depending on the operating region 
of the system. 

In conclusion, the value of a good h will depend on the channel codes, the relay position 
and channel conditions. This corroborates the previous statement that a random choice of the 
NC-coefficient might lead to some performance degradation. 

Finally, note that if the channel codes are composed of several subgraphs (e.g. LDPC or 
Turbo codes), several transfer functions (one for each compounding subgraph) are obtained 
and a direct representation of these functions will result in an A^-dimensional EXIT chart. By 
converting the A^-dimensional EXIT chart into a two-dimensional EXIT chart, for example 



by using the EXIT Chart Projection Algorithm proposed in p9| |, the above analysis could 
still be applied. 

V. Simulation Results 

In order to asses the performance of the proposed scheme and to corroborate the conclu- 
sions from the previous analysis, several sets of simulations have been performed. We have 
considered a symmetric scenario where both source nodes are placed at the same distance 
from the destination nodes, i.e. diD = d2,D- Regarding the relay position, three different 
scenarios are defined: 



The minimization is performed by an exhaustive search over all possible h G GF{2'') (Figure 6 
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A) The relay and the sources are deployed at identical distance from the destination, (i.e., 
SNRr^d = SNRi^D + OdB). This setup was used in 0, Q. 

B) The distance between the relay and the destination is three quarters of the source- 



destination distance (SNRr^d = SNRi^d + 4.4dB) as used in [[T0|, [[17 1. 
C) The relay is set at approximately half the distance between the sources and the destination 

(SNRr^d = SNRi^o + lOdB). This setup was also considered in Q. 
The channel codes used at both source nodes are identical, non-systematic convolutional 
codes of rate 1/3. Two types of convolutional codes are considered: a 2^-state [5, 7, 7]8 codqj 
(heretofore denoted as CC2), and a 2^-state [554, 624, 764]8 code (correspondingly, CC6). A 
zero-bit termination tail is appended at each source sequence. We use packets of i^ = 32 
bits (i.e. we use M = 32 ■ 3/2 = 48), and the interleavers have been randomly generated 
and are independent of each other. We have considered 4-QAlVJj which leads to an spectral 
efficiency of p = 4/9 [bits per complex dimension], and the number of iteration for the 
SPA has been set to X = 15. Finally, from the region of achievable decode-and-forward 



rates of the TD-MARC given in [18|, (derived from the capacity bounds of the C-MARC in 
||3|), a set of upper bounds on the outage probabilities has been obtained by specifying the 
actual packet lengths of the source information bits and the transmitted sequences (K = 32, 
M = 48). These outage probability bounds are used as an information-theoretic benchmark 
for the Packet Error Rate (PER), where a packet is in error if one or both sources packets 
are erroneously decoded. 



Based on the analysis carried out in Section IV the NC-coefficients for g = 3 that 



minimizes P{T^(1) < x} and V{T^{0) < x} are given by h = (6,6) and h = (1,1), 



The subindex 8 in the definition of the code stands for octal 

^Further improvement can be expected with higher order modul 
Coded Modulation (BICM |30|) technique. However, this research line lies beyond the scope of this manuscript. 



^Further improvement can be expected with higher order modulation schemes and by applying the so-called Bit-Interleaved 
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respectively. To corroborate the optimality of these values a first set of simulations have been 
done for scenarios A and C (Fig. |7]j8]). Since we are mainly interested in the selection of the 
NC coefficients, in these simulations error-free links between the sources and the relay are 
assumed. Nevertheless, the noisy source-relay link case is also discussed later in this section. 
By using the channel code CC2, Figure [v] plots the end-to-end PER versus Eb/No = 
SNR — lOlogj^oP (i^ dB) for all possible values of h and for scenarios A (lefthand plot) 
and C (righthand plot). It can be observed that in scenario A, and for all values of SNR, the 
coefficients h = (6, 6) (i.e. h^ are the optimal choice. For scenario C, h^ is still optimal 
in the range SNR > SNRth (recall that SNRth denotes the SNR at the crossing of the 
curves for h^ and hg); however, for SNR < SNRth the optimal choice is h = (1,1) (i.e. 



hg). This corroborates what was stated in Section IV that the use of low-complexity codes 
reduces the influence of Pjearly-crossing}, leading to low values of SNRth (in scenario A 
SNRth < and in scenario C SNRth ~ 3). Moreover, it can be seen that as the quality 



of the relay-destination link increases (going from scenario A to C), the value of SNR 



th 



increases; as also stated in Section IV 



Figure [8] shows a similar analysis on the PER when using the CC6 code, instead. It can 
be observed that in both scenarios SNRth arise at high values of SNR, since hg outperforms 
h^ for the SNR ranges of interest (in scenario A and C, the crossing occurs at SNR > 20dB 
and SNR > 16dB, respectively), corroborating the fact that we are using high-complexity 



convolutional codes (see Section IV). 



Furthermore, these results show that as the relay is placed toward the destination, a tailored 
selection of the NC-coefficients improves the performance of the system; e.g. in scenario 
C, gains of 1.3dB and 1.5dB, with respect to the worst selection, can be achieved when 
using CC2 and CC6, respectively . One can also conclude from the figures that, since the 
performances for the different coefficients are roughly uniformly distributed between the best 
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Scenfuio A; SNBjw = SNR_^!j + f)dB 



Scenario C^ SNHro = SNRsd + lOdD 




Fig. 7. Outage probabilities and PER performance of the proposed approach when the CC2 convolutional code is used. 



Scenario C: SNRrd = ■'^NRsD + lOrfS 




Fig. 8. Outage probabilities and PER performance of the proposed approach when the CC6 convolutional code is used in 
scenario A (left) and C (right). 



and worst case scenarios, a gain of 0.75dB (scenario C and CC6) is obtained with respect to 
the random choice of coefficients. 
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The next set of simulation compares the gap between the performance of the proposed 
system with respect to the outage probability. This gap is also computed for different schemes 
found in the literature, all having a spectral efficiency less than 2 (bits per complex di- 
mension). It should be mentioned that all these schemes [[6), [[8j, [10|, [17| were analyzed 
assuming error free source-relay channels. Therefore, the outage probabilities (upper bounds) 
for scenarios A,B,C has been computed with the later assumption and for unconstrained 
channel inputs (Gaussian). 

These gaps are given in Table |I} where the first two columns refer to our scheme. Observe 
that the proposed scheme outperforms all the reference schemes. Moreover, the proposed 
system achieves these results by using short-length codewords: 144 complex dimensions per 
use of TD-DF-MARC in contrast to the 2176 utilized in [[T0|, the 6000 utilized in Q, [|8| 



or the 27000 utilized in [17| 



TABLE I 
Gaps to the outage probabilities 



Scenarios 


q=3, CC2 


q=l, CC6 


[ 


17 


\ 


[ 


6 


\ 


[ 


8 


1 


(10 


1 


Scenario A 


1.39 dB 


1.44 dB 


- 


3.4 dB 


2.7 dB 


- 


Scenario B 


1.64 dB 


1.64 dB 


1.7 dB 


- 


- 


4.8 dB 


Scenario C 


2.04 dB 


2.36 dB 


- 


- 


5.2 dB 


- 



Finally, we consider the case where non-ideal source-relay channels are used. An un- 
changed implementation of the proposed scheme will lead to error propagation at the decoder 
and as a result a degradation on the performance as shown in Figure |9} Moreover, at high 
SNRs, the closer the relay to the destination is, the stronger the impact of the error propagation 
will be. However, since the error propagation also degrades the outage probability as shown 
in Fig. |9} the gap between the performance of the unchanged scheme and the new outage 



probability curve [18| is still small, and is in fact negligible at high SNRs. We conclude that 
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our scheme is still robust when non ideal channel error propagation occurs. 

An alternative scheme to cope with the error propagation at the destination is as follows. 
This scheme assumes that the destination node knows if an error has occurred at the relay 
(Relay State Information, RSI), e.g, by using an additional low rate error free relay-destination 
channel. To this end, a cyclic redundancy check (CRC) is added at the end of the source 
transmitted sequences, so that the relay can detect if any residual errors have occurred in X^ 
and X^. If errors are detected in X™, the relay sets /i™ to (i.e. the sequence transmitted 
by the relay only conveys information from 5^~*) and the relay communicates that an error 
has been produced to the destination. The decoder sets at the network check node, h^ to 0, 
avoiding in this way error propagation through channel decoders. However, some performance 
degradation is still expected since p^ARc^^m^ ^^^^ p^g ^ j^ ^^^ uniformly distributed and 
therefore no iteration gain will be obtained. The performance of the RSI scheme is shown 
in Figure |9] together with the corresponding outage probability. One can observe that the 
performance loss with respect to the error-free case is negligible. 

VI. Conclusions 

This paper proposes a novel joint network-channel coding scheme for the Decode-and- 
Forward Time-Division Multiple Access Relay Channel. Specifically, we have designed a 
Joint-Network-Channel code which does not perform channel coding on the already network- 
coded bits, reducing the complexity at the relay node without compromising performance. 
A method for selecting the best pair of coefficients of the linear combination is derived 
based on an EXIT charts analysis. Moreover, the proposed code allows the sources to use 
completely different channel codes, at the sole expense of an increased complexity on the 
EXIT chart analysis. The decoding at the destination is performed by applying the SR\ 
over the derived factor graph of the JNCC code. Monte Carlo simulations show that the 
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Fig. 9. PER performance of the proposed schemes for the erroneous source-relay channel setup when the CC2 convolutional 
code and h — (6, 6) is used. 



proposed scheme outperforms, in terms of its gap to the corresponding outage probability 
upper bounds, the previously published schemes for the same network setup. Besides, this 
gain is achieved by using short-length codewords, which makes our proposal particularly 
attractive for low-latency applications. 
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